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Abstract 

We introduce an extension of the diagrammatic rules in random matrix theory and 
apply it to nonhermitean random matrix models using the approximation. A 
number of one- and two-point functions are evaluated on their holomorphic and non- 
holomorphic supports to leading order in 1/N. The one-point functions describe the 
distribution of eigenvalues, while the two-point functions characterize their macro- 
scopic correlations. Generic form for the two-point functions are obtained, gener- 
alizing the concept of macroscopic universality to nonhermitean random matrices. 
We show that the holomorphic and nonholomorphic one- and two-point functions 
condition the behavior of pertinent partition functions to order 0{1/N). We derive 
explicit conditions for the location and distribution of their singularities. Most of 
our analytical results are found to be in good agreement with numerical calculations 
using large ensembles of complex matrices. 
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1 Introduction 



Random matrix models provide an interesting setup for modeling a num- 
ber of physical phenomena, where noise plays a prominent role. Applications 
range from atomic physics to quantum gravity. In these phenomena, one may 
distinguish between the universal behavior of noise exhibited as fluctuations 
at the mesoscopic scale (microscopic universality) and the macroscopic scale 
(macroscopic universality) in bulk systems with or without absorption. Sim- 
ple physical realizations of these universalities are encountered in disordered 
mesoscopic systems [1]. 

In this work, we will be concerned with the generic behavior of two classes of 
random matrix models: those without dissipation (hermitean) and those with 
dissipation (nonhermitean) . Since conventional hermitean matrix ensembles 
have been considerably investigated [2] (For recent discussion on macroscopic 
correlators see e.g. [3]), we will keep our discussion to the chiral unitary en- 
sembles [4], and we devote most of our studies to nonhermitean ensembles, 
which are less known although they play an increasingly important role in 
quantum problems [5,6]. 

A number of methods have been devised to calculate with random matrix 
models. Most prominent perhaps are the Schwinger- Dyson approach [7] and 
the supersymmetric method [8]. In this paper, we would like to elaborate on 
diagrammatic techniques for a variety of random matrix models, with special 
focus on the chiral and nonhermitean ensembles. The use of diagrammatic 
methods in random matrix models is not new, although those methods were 
not previously applied to nonhermitean matrices, and we refer the reader to 
[9] (and references therein). In this work we will discuss one- and two-point 
functions. They carry information on the eigenvalue distributions and their 
(macroscopic) correlations. They involve singular operators in the large N 
limit, and in general break spontaneously holomorphic symmetry. We show 
also that they condition in a simple way the character of pertinent partition 
functions to order 0{1/N), as in general expected. 

The organization of this paper is as follows: in section 2, we set up the notations 
for the hermitean case, and evaluate the resolvent in the approximation. 
We show how the result changes in the chiral case, and apply it to a schematic 
version of a Dirac operator. In section 3, we extend our diagrammatic argu- 
ments to the nonhermitean case, using real and nonhermitean random ma- 
trix models as an example. The important difference between hermitean and 
nonhermitean matrix models, is the possibility of spontaneous breaking of 
holomorphic symmetry in the eigenvalue plane. In section 4, we work out the 
resolvents in the case of real and asymmetric matrices, nonhermitean chiral 
matrices and a random scattering model. In section 5, we evaluate a number of 
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two-point correlators. We show that they all follow the general lore of macro- 
scopic universality, thereby generalizing the idea of Ambj0rn, Jurkiewicz and 
Makeenko [10] and also Brezin and Zee [11] to the nonhermitean case. All our 
results are found in good agreement with large scale numerical calculations. 
In section 6, we show that the one- and two-point functions may be used to 
generate pertinent partition functions to order 0{1/N) for both the holomor- 
phic as well as nonholomorphic case. We use the Ginibre-Girko ensembles [12], 
as a way of illustration. A Straightforward generalization to other unitary en- 
sembles is briefly mentioned. Our conclusions are summarized in section 7. 



2 Hermitean Random Ensembles 

To illustrate the arguments and set up the notation for the diagrammatic ar- 
guments to follow, let us first consider the well known case of a deterministic 
plus random hermitean ensemble with Gaussian distribution. It is convenient 
to use the diagrammatic notation introduced by [11], borrowing on the stan- 
dard large N diagrammatics for QCD [14]. For the random part, the fermionic 
Lagrangian is simply 

C^U<-H'M\ (2.1) 

where if is a hermitean random matrix with Gaussian weight. We will refer to 
■0 as a "quark" and to if as a "gluon" . The "Feynman graphs" following from 
(2.1) allow only for the flow of "color" (no momentum), since (2.1) defines a 
field-theory in + dimensions. They are an efficient way of keeping track of 
"color" factors. ^ The "quark" and "gluon" propagators (double line notation) 
are shown in Fig. 2.1. For hermitean random matrices the Gaussian averaging 
preserves holomorphic symmetry in the large N limit. This is not the case for 
nonhermitean random matrices when singular operators are considered (see 
sections 3-6). 




1 

z 



Fig. 2.1. Large N "Feynman" rules for "quark" and "gluon" propagators. 



The names "quarks", "gluons", "color" etc. are used here in a figurative sense, 
without any connection to QCD. In order to avoid any confusion, we always use 
these names in quotation marks. 
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2.1 Pastur Equation 



The pertinent resolvent for (2.1) is given by [11,15,16] 



where Ti, is an N x N hermitean random matrix, and D a deterministic part. 
The distribution of eigenvalues ^{z) oi H — D + Ti follows from (2.2) through 
7ru{z) — dzG. 

Throughout, the weight average is considered Gaussian for simplicity. How- 
ever, all our arguments generalize to the non-Gaussian weight. In terms of the 
one-particle irreducible self energy E, (2.2) reads 

G{z) = \-tY ^ ^. (2.3) 

N z-D-H ^ ' 



We note that since the argument in (2.2) displays poles on the real axis only 

(hermitean matrices), it is safe to assume that after averaging and in large A^, 
the resolvent is holomorphic in the 2;-plane modulo possible singularities on 
the real axis, for which v{z^ 7^ 0. This property is manifest in (2.3). 



Fig. 2.2. Rainbow diagrams. 

In the large limit the equation for the self energy E follows from resum- 
ming the rainbow diagrams of Fig. 2.2. All other diagrams (non-planar) are 
subleading in the large N hmit. The result is 



where T>ah,cd is the random matrix propagator given by 

2^a6,c<i =< nlnl >= j^SadSbc ■ (2.5) 



Inserting this into (2.4) we find that E = E • 1, where 1 is the identity matrix, 
and E is a scalar satisfying Pastur's equation (see Fig. 2.3) 

— tr = E. (2.6) 

N z-D-E ^ ' 
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We note that S = ^i^)- 1^ is convenient, for the generaUzation to the chiral 
case, to represent the propagator as a tensor product of matrices 



P= 1® 1 



(2.7) 



As a result equation (2.4) becomes 



1 

TV" 



-tr 



1 



z-d-j: 



1 = E 



(2. 



This result was derived by a number of authors [11,15,16]. 




Fig. 2.3. Pastur equation. 
2.2 Chiral Pastur Equations 

In the chiral case the + dimensional Lagrangian (2.1) becomes 

£ = Mz5', - MDi^' (2.9) 



with the chiral random matrix 



M 



n 
7^t 



(2.10) 



Here, the "quarks" are chiral in the sense that ip = {iPr,iPl), following the 
block decomposition in (2.10). The chirality operator is 



75 



1 
-1 



(2.11) 



We may now apply the results of the previous section but with a different 
propagator. Namely 



T>ab,cd =< M^M'a >= -5ad5bc a> N,b<N or a<N,b> N (2.12) 
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and zero otherwise. In the tensor product notation we get 




(2.13) 



Again, the random averaging is carried using a Gaussian weight. Since (2.9) is 
chiraUy symmetric for 2; = 0, the question arises as to whether the averaging 
breaks spontaneously chiral symmetry in the large N limit. This is indeed the 
case for 2; — > iO after — > 00. Inserting (2.13) into equation (2.4) we obtain 





(2.14) 



We see that S is diagonal 



^ 



Eo- 1 



(2.15) 



with El and E2 scalars and the Green's function is G{z) — |(Ei + E2). The 
chiral Pastur's equations now take the form 



1 



1 



z-D 

1 




Si, 



= Eo. 



(2.16) 
(2.17) 



As z — > iO, El = E2 7^ implying that in the large N limit, the Gaussian av- 
eraging breaks spontaneously chiral symmetry. This is particularly interesting 
as it shows the following: the mechanism of spontaneous symmetry breakdown 
requires minimally that the interaction is chirally symmetric and random, a 
point presently embodied in the instanton vacuum description [17]. A related 
idea was first discussed in [18]. 



7 



2.3 Random Dirac Operators 



As an illustration of the above equations, consider the schematic case of chiral 
fermions with mass m. In this case, the deterministic 2N x 2N matrix is 



D = 




(2.18) 



Rearranging the rows and columns of the matrix z — D — f] and inverting, 
yield the following equations for Ei and E2 



z-m-Ei 



(z — m — Tji){z + m — E2) ' 

=^ = 7 'trf yy (2.19) 

[z — m — ^i)[z + m — 

Using the fact that E2/E1 = {z + m)/{z — m) and G — (Ei + E2)/2 we get 
the equation for the Green's function 

{G^ - 2zG^){m^ - z") + Gz^{m^ - {z" + 1)) + / = . (2.20) 



We note that when m = this gives exactly z^ times the ordinary - hermitean 
- Pastur's equation for this case. Solving for Ei gives 

G=-fl + ^5^^1 (2.21) 
2 y \J z^ — I 



which is holomorphic in the 2;-plane modulo cuts on the real axis. The solutions 
to (2.19) were used in [19]. 



3 Nonhermitean Random Ensembles 



If we were to use nonhermitean matrices in the resolvent (2.2), then configu- 
ration by configuration, the resolvent displays poles that are scattered around 
z = Q (for L> = 0) in the complex 2;-plane. Their spread is given by the 
Gaussian distribution, which is of order one. In the large limit, the poles 
accumulate in general on finite surfaces (for unitary matrices on circles), over 
which the resolvent is no longer holomorphic. The (spontaneous) breaking of 
holomorphic symmetry follows from the large N limit. As a result dG/dz 7^ 
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on the nonholomorphic surface, with a finite eigenvahic distribution. In this 
section we wiU set up the diagrammatic rules for investigating nonhermitean 
random matrix models. 

For hermitean matrices, "quarks" and "conjugate-quarks" decouple in the 
"thermodynamical" limit {N oo). Their respective resolvents follow from 
the + dimensional Lagrangian 

Co^^{z-M)7j; + ^{z-M^)(l) (3.1) 



and do not 'talk' to each other. They are holomorphic (antiholomorphic) func- 
tions modulo cuts on the real axis. For nonhermitean matrices, this is not the 
case in the large N limit. The spontaneous breaking of holomorphic symmetry 
in the large N limit, may be probed in the 2;-plane by adding to (3.1) 

Cb = A(^0 + 0^) (3.2) 



in the limit A — > 0. The combination Cq + Cb will be used below to in- 
vestigate a number of random and nonhermitean ensembles with or without 
deterministic parts. 

From (3.1,3.2) we define the matrix- valued resolvent through 



\Qqq Qqq 



z-M 
A 




(3.3) 



where the limit A?" — > oo is understood before A — 0. For convenience, the 
"quark" field is organized in an isodoublet ^ = ("0)0) — (?)?)-The "quark" 
spectral density follows from Gauss law, 

iy{z, z) = -d, G{z, z) = -\-d, Tr^ G,, (3.4) 



which is the distribution of eigenvalues of Ai. For hermitean Ai, (3.4) is val- 
ued on the real axis. As A ^ 0, the block-structure decouples, and we are left 
with the original resolvent. For z — > +iO, the latter is just a measurement of 
the real eigenvalue distribution. For nonhermitean (3.4) is valued in the 
2;-plane. As A ^ 0, the block structure does not decouple, leading to a non- 
holomorphic resolvent for certain parts of the 2;-plane. Holomorphic symmetry 
is spontaneously broken in the large A^ limit. The size of the nonholomorphic 
regions is conditioned by the divergence of the "quark-conjugate-quark" two- 
point function (see below). 
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3. 1 Example 



To set up the diagrammatic rules we will use an example. Consider the case 
of real and nonhermitean matrices with off-diagonal correlations. Specifically, 
let M be a real and asymmetric N x N matrix so that 

(Mj) = 1 {M,jM,,) = ^ . (3.5) 

The case r = 1 reduces to symmetric real matrices. The random ensemble 
described by (3.5) occur naturally in neural network problems, and have been 
discussed in [20] using the replica construction. 



3.2 Diagrammatic rules 



From (3.1,3.2) wc note that there arc two kinds of "quark" propagators {1/z 
for "quarks" ip and 1/^ for "conjugate-quarks" 0, where both can be incoming 
and outgoing) and four kinds of "gluon" propagators, associated to the fol- 
lowing "quark-conjugate-quark" amplitudes (treating the "quarks" as external 
sources) 



(0,(M^)^0>,(M^)S0'') =0,0>,0''((M^)^(M^)S) . 
In each of these four cases, averaging is done with the weight 



(3.6) 



<...>= 



j[dM\ 



exp 



N 



2(1 -r^) 

reproducing the correlations (3.5). Indeed, 



tr(MM^ - tMM) 



(3.7) 



< M^M^ >-- 



1 

TV 



(M'^'^ + r^.^^c") 



(3.^ 



where the second contribution corresponds to twisting the lines with a "penalty 

factor" r. For r = ±1 this is equivalent to symmetrizing or antisymmetrizing 
the double lines in the "gluon" propagator. All combinations arising from the 
averaging (3.7) are shown in Fig. 3.1. The pairs (ae),(bf), (eg) and (dh) in 
Fig. 3.1 correspond to four "gluon" propagators, evaluated with the weight 
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(3.7). The lower row in Fig. 3.1 corresponds to terms with r 7^ 0. The arrows 
denote incoming and outgoing particles, the "quarks" and "conjugate-quarks" 
are distinguished by their labels. 



a) 



q q 
q q 

X 

q q 



■a 
■d 



b) 




q q 



c) 




q q 



^ ^ ^ ^ 



d) 



q q 



C— ^ 



q q 



, q q 

c^*-l^^d 

q q 



Fig. 3.1. All "gluonic" amplitudes for real nonhermitean matrices. 

The generalization to the hcrmitcan-antihcrmitcan case, i.e. the case of com- 
plex matrices, corresponds to the substitution with the measure 



<...>= j [dM] exp 



(1 - r2) 

The relations (3.5) should be replaced by 



tr(MMl' - tRqMM) 



(3.9) 



1 

TV 



N 



(3.10) 



corresponding to hermitean {t — 1), antihermitean (r = —1) or general com- 
plex (r = 0) matrix theory. In the latter the relevant "gluonic" amplitudes 
correspond now to Fig. 3.1a-3.1d. 



4 One-Point Functions 



In this section we will explicitly calculate a number of one-point functions 
(resolvents) for several nonhermitean random matrix models, using and ex- 
tending when needed the diagrammatic approach discussed above. 



Pure Complex Random Matrices 



We seek the one-point functions for (^j^j) as a 2 x 2 matrix. The equation for 
the one particle irreducible (IPI) self-energy follows from Figs. 3.1-4.1 in the 
form 
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^ El E2 ^ 

y S3 S4 y 



N 



tr. 





(r 








;) 




.1 





El 


A 




-1 















S3 


z 


— E4 ^ 




[ir) 



(4.1) 



Here the trace is meant component-wise (block per block), and the argument 
of the trace is the dressed propagator. The operation o is not a matrix multi- 
plication, but a simple multiplication between the entries in the corresponding 
positions. Here tr^r is short for the trace on the N x N block-matrices. 

The first correlator in (3.5) does not influence the "quark-quark" interaction 
- it corresponds to the double line with a twist in the corresponding Pastur 
equation - hence subleading. However, this twist could be compensated by 
the twisted part of the propagator coming from the second correlator (3.5), 
thereby explaining the factor r in the upper left corner of (4.1). This is shown 
in Fig. 4.1, where the first (subleading) graph comes from Fig. 3.1a and the 
second (leading) comes from Fig. 3.1e. The other entries in (4.1) follow from 
Fig. 3.1 by inspection. 



The "quark" one-point function is now 

1 



N 



tfN Sqq = (^ - ^4)/det. 



(4.2) 



It follows that E2 = E3, with 



det El = r{z — E4) 
det E4 = t{z — El) 
det E2 = E2 — A , 



(4.3) 
(4.4) 
(4.5) 



where det = {z ~ ^i){z — E4) — (A — E2)^. Substituting r = E2 — A in the last 
relation in (4.5) yields the equation 



((^-Ei)(5-E4)-r')(r + A)=r. 



(4.6) 



For A = 0, the solution with r = is holomorphic while that with r 7^ is 
nonholomorphic. In the holomorphic case, Ei(2; — Ei) = r, and the resolvent 
is simply 



G{z) 



2^ 



(4.7) 
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where the upper sign corresponds to the solution with the pertinent asymp- 
totics. In the nonholomorphic case, G{z, z) — z — E4, with 



G{z,z) 



Z — TZ 



(4. 



in agreement with [20] . The boundary between the holomorphic and nonholo- 
morphic solution follows from the condition \G{z,z)\'^ — \G{z)\'^ — 1, that 
is 



X 



;1 + t)2 {l-Tf 



1 



(4.9) 



which is an elhpse. Inside (4.9) the sohition is nonholomorphic and outside it 
is holomorphic. The case investigated by Ginibre and Girko [12] follows for 
T = 0. The one-point function for the random matrix [Hi -\- iH2)/\/^ can be 
obtained similarly. The self-energy equation being the same as before but with 

T = 0. 




q q ^ q q q q ^ q q 

Fig. 4.1. Self-energy equation for real nonhermitean matrices. 



4-2 Real Asymmetric (Square) Matrices 

Following Khoruzhenko [21] we define the random matrix ensemble 

Ho + A (4.10) 



where Hq is a deterministic but otherwise arbitrary (e.g. complex) matrix and 
A is a real asymmetric random matrix with 



(4.11) 



The IPI self energy equations are now 







= tTN j 


\ S3 







z — Hq — A — 



A-E, 











^0 


vA 















(4.12) 
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Care must be taken with the definition of the inverse matrix. The inverse is 
taken only with respect to the "isospin" indices (2x2), so the determinant is 
now a matrix in N-space. We obtain 





( 


Si J 





- Hq — XI4 S2 — A 

E3 — A z — Hq — Eiy 



( v'^D-/ 



v^Dl' 



(4.13) 



with E2 = E3 = X for A and 



DL = iz- Hi - ^,){z -Ho- El) - , 

Dr^{z-Ho- Ei)(z - Hi - E4) - (4.14) 

since we do not impose that Hq and hI commute. We see that Ei = E4 = 
and the equation for E2 = x now has the form 

x^v'-x-j^ tTMiiz - Hl){z - Ho) - xY' . (4.15) 



The solution with x 7^ corresponds to the non-analytical Green's function. 
Its "quark" Green's function is given by 



G{z,z) = ^tr^[{z-Hl)Go{-x') 



(4.16) 



where Go{w) is the generalized resolvent 



Goiw) = 



{z-HI)(z-Ho)+w 



(4.17) 



and X satisfies 



(4.18) 



These are the equations obtained by Khoruzhenko [21]. 



4.3 Nonhermitean Chiral Matrices 



In this section we will reanalyze the chiral random matrix model of section 2.2 
in the presence of a nonhermitean part, e.g. "chemical potential" jj,. The moti- 
vation for that stems from the constant mode sector of the massless and chiral 
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Dirac operator at finite chemical potential [22]. Here we will use Gaussian 
chiral and random matrices as described by (3.1). 



If we were to define 



7 = ^7o 



-1 

1 



and M 



' -n + n 
^n + n'' 



(4.19) 



then the IPI self-energy equations in the planar approximation are given by 



El E2 
y S3 S4 



1 I z - ii'j -El A - E2 
Tftrjv 

' A - E3 z + jd^-E^ 



(4.20) 



where T> is the "gluon" propagator (2.12), and Ej are diagonal 2N x 2A'" 
matrices. Inverting in (4.20) with respect to the "isospin" indices gives 







= tr^v 






V 



2; + /i7 - E4 E2 - A 
\^ E3 - A z - /i7 - El 



V V 

V V 



■ (4.21) 



Here A is the determinant of the 2x2 matrix with 2N x 2A'" entries. The 
equation for E2 7^ in (4.21) gives 



l = ltr (A-.P) 



(4.22) 



with 



A-i 



-^{z — E4 — z + Ei)/q; 



^IJi{z - E4 - ^ + Ei)/a 



(4.23) 



and where a^^ is the determinant of the 2A'" x 2N matrix A. Inserting (4.23) 
into (4.21) yields 



E4 = 2; — El — /i 
El = z — E4 + 



2Z — E4 — z -\- El 



a 

2^ — E4 — 2; + El 
a 



(4.24) 
(4.25) 
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In particular Ei + E4 = Re^; = x. The resolvent for the "quark-quark" part 
is simply Ei. Solving (4.25) in conjunction with (4.22) gives 

q;2-q; + //^(^-E4-z + Ei)2 = 0. (4.26) 



Simple algebra for the "quark-quark" resolvent gives 



Gfe^-) = ^tr„e„ = £-i,-i-^ (4.27) 



a result first derived in [23] using different arguments. 

For E2 = we recover the holomorphic solution [23,34], Ei(2;) = G{z)l, 
E4 = E|, with G{z) fulfilling the cubic Pastur equation 

G%z) - 2zG\z) + {z^ + + l)G{z) -z^O. (4.28) 



4-4 Random Scattering Model 



As a final application of the techniques introduced above, we consider the 
random matrix model discussed by Mahaux and Weidenmiiller [24] in the 
context of chaotic resonance scattering. The model is generically described by 
a nonhermitean Hamiltonian of the form 

H - i-fVV^ (4.29) 



where H is random Gaussian (orthogonal) and N x N matrix- valued, while V 
is N X M matrix- valued [6,25]. The form (4.29) is dictated by unitarity. The 
V propagator is given by 

(VCVt) = Is'^Xi . (4.30) 



ak I b 

q q 



ak I b 

q q 



a,b=l,...,N 
k,l = l,...,M 



Fig. 4.2. Vertex corresponding to the VV^ terms in the Hamiltonian. 

In the following we will use the notation 6 = i'~f and m = M/N. The quadratic 
appearance of V in (4.30) calls for additional changes in the preceding Feyn- 
man rules. In particular, the IPI self-energy has two new contributions: 
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(1) a channel- vertex correction: 

dm for "quark" , —5m for "conjugate-quark" , and for mixed. 

(2) a sum of rainbow diagrams 

each of which is shown in Fig. 4.3. Here the extra diagrammatic comphcations 
stem from the fact that two V "gluon" hues couple to either the "quark" or 
the "conjugate-quark" lines through the new vertex (see Fig. 4.2). 



, " 

mFqq 




1^1 q q 
q q 

S2=/^ + 

Q Q q q 

Fig. 4.3. Self-energy equations for the random scattering model. 
The self-energy equations are explicitly given by {G — Gqq, G — Ggg) 



El = G + mFqq , 
S4 = G + mFqq , 

^2=Gqq+mFqq, (4.31) 
S3 = Gqq + mFqq 

where G — {z — T,4)/det, Gqq — (E2 — \)/det, etc., with det — {z — Si)(^ — 
S4) — (A — S2)^. The equations for the F's are just the generating ( "Lippmann- 
Schwinger") equations for the rainbow graphs, that is (see Fig. 4.4) 



Fqq — 6 + SGFqq + SGqgFqq , 
Fqq — —5 — 6GFqq — 6GqqFqq , 

Fqq = -SGFqq - SGqqFqq , (4.32) 
Fqq — SGqqFqq + SGFqq . 
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One may solve for Fqq to obtain 



denom 



The equation for the mixed self-energy E2 7^ now gives 



(4.33) 



where denom is 



det^l- (4.34) 
denom 



- 5^ 

denom = 1 + 6{G-G) - —. (4.35) 

det 



q q qq q qq q <i q qq q <i 




7 

qq q qq q ^ q qq q ^ 

q qq q q q qq q ^ 

Fig. 4.4. "Lippmann-Schwinger" equations for F's. 

The way to proceed now is as follows. One can write an equivalent generating 
equations for Fqq by adding the graphs on the "quark" side. This shows that 
Fqq — Fqq. This iu conjuuctiou with (4.34) gives the relation 

Fm + Fqq = -{G + G){1 - det) . (4.36) 
m 



Inserting this into the self-energy equations one obtains 

G + G = x (4.37) 



Re-expressing Fgg by Fqq in the appropriate self-energy equations and using 
(4.34) one obtains 

z-Gdet + 5zG^5(l-m) + G{l + 6x) - 6 /det . (4.38) 
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Eliminating G, one solves this in conjunction with (4.34) treated as an equa- 
tion for the determinant det: 

[1 + 5(x - 2G) - 5"^ /det] det ^ 1 + 5(x - 2G) - 5'^/det - . (4.39) 



The final result (after changing 5 into — ^7) is 



7 



1 m 

7 y ^-ly 



(4.40) 



This result was first derived by Haake, Sommers and coworkers using the 
replica method [25], and the supersymmetric method [26]. 

Before closing this section we would like to mention briefiy the connection of 
our diagrammatical approach to the mathematical concepts of free random 
variables [13,27]. For applications of free random variables to various physical 
problems see e.g. [28-32]. With the compact notations 





Y Gqq Gqq 



the equations (4.31) and (4.32) could be rewritten in the form 



(4.41) 



(4.42) 



Solving for the matrix S we get 



(4.43) 



Introducing the generalized Blue's function [33] (an extension of Zee's ap- 
proach [31] to nonhermitean ensembles), as 2 x 2 matrix valued function de- 
fined by 



B{G) = Z 




{AAA) 
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with A ^ 0, we see that S = B(^) - g-^. Therefore 



Big)^g + g-' + mii-5g)-^5 (4.45) 

in which we recognize the generahzed addition law postulated in [33], with 
generalized Gaussian Blue's function 



Bgauss{A)^j + A (4.46) 

and generalized Blue's function for the random scattering part of the Hamil- 
tonian form 

B_i^vv^ (A) = m(l - SA)-^S + (4.47) 



with 




(4.48) 



This completes the diagrammatic proof of the result announced already in our 
earlier work [33]. 



5 Two-Point Functions 



To probe the character of the correlations between the eigenvalues of non- 
hcrmitean random matrices, cither on their holomorphic or nonholomorphic 
supports, it is relevant to investigate two-point functions. A measure of the 
breaking of holomorphic symmetry in the eigenvalue distribution is given by 
the connected two-point function or correlator 

)c (5.1) 



tr 



1 

z-D-n 



where the z and z content of the averaging is probed simultaneously. In [34] 
the correlation function (5.1) was shown to diverge precisely on the nonholo- 
morphic support of the eigenvalue distribution, indicating an accumulation in 
the eigenvalue density. In the conventional language of "quarks" and "gluons" , 
(5.1) is just the correlation function between "quarks" and their "conjugates". 
A divergence in (5.1) in the z-plane reflects large fluctuations between the 
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eigenvalues of the nonhermitean operators on finite z-supports, fience tlieir 
"condensation" . 

It was sfiown in [10] and [35] tfiat for hermitean matrices (with z — > ty) the 
fluctuations in connected and smoothed two-point functions satisfy the gen- 
eral lore of macroscopic universality. This means that all smoothed correlation 
functions are universal and could be classified by the support of the spectral 
densities, independently of the specifics of the random ensemble and genera 
in the topological expansion (see [36] for a recent discussion). In this section 
we will evaluate a number of two-point functions for random but nonher- 
mitean ensembles on both their holomorphic and nonholomorphic supports. 
We confirm that the general lore of macroscopic universality extends to the 
nonhermitean case as well. 

The concept of smoothed correlation functions will become manifest when 
comparing to numerical calculations. Indeed, when correlating N xN matrices 
at finite eigenvalue-separations, oscillations are expected due to the occurrence 
of a large number of eigenvalues. The correlation functions produced in the 
1/N analysis are smoothed. Unsmoothed correlation functions are of interest 
for the studies of spectral form factors in the crossover region [37] . They will 
not be discussed here. 



5. 1 Notations 



To establish the notations for the two-point correlators, we follow [35] and 
define the two-point connected (c) correlator for GUE, 



^ / x / 1 11 1 
GJz, w) — ( — tr tr 

^ ' ' \n z-HN w-H, c 

= 9A(^trlog(z-?i)ltrlog(«;-7i)) . (5.2) 



c 



Expansion of the logarithms yields 



Let us first consider n = m. The graphs in this case are represented by the 
wheel diagram (Fig. 5.1a), where the rungs of the wheel are given by the 
"gluon" propagator from Fig. 2.1. The diagram splits into n disconnected 
sectors, corresponding to the fact that < trTi^trH" >— n. Resumming the 
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diagonal terms give to 

N''Gc{z,w) = 



(5.4) 



All the remaining (planar) diagrams are rainbow-like (see Fig. 2.2) and corre- 
spond to dressing the bare "quark" propagators, i.e. l/z — > G{z). The result 
is 

N'^Gciz, w) = -d,d^ log[l - G{z)G{w)\ . (5.5) 

The two-point correlator follows functionally from the one-point correlator. 
Such a relation for the case of Gaussian ensemble was already noted in [38] 
prior to the universality arguments. 

This result can be further reduced [10,35], to show that the two-point cor- 
relator depends only on the end-points of the spectral density irrespective 
of the choice of the weight. Since the smoothed correlation function between 
the eigenvalues follow from (5.5), it reflects on correlations between eigenval- 
ues separated by in the spectrum. Its generic form reflects on macroscopic 
universality in the eigenvalue correlations. This is to be contrasted with mi- 
croscopic universality which is a statement about the correlations between 
eigenvalues one- level spacing apart in the spectrum [2] . 




Fig. 5.1. Sample wheel diagram (n = 4) for the two-point correlators (GUE), (a) 
and generalized wheel diagram with kernel V (b). 

5.2 Ginibre-Girko Correlator 

In this section we show how to generalize (5.5) to the nonholomorphic region 
of complex and random matrices. We begin our discussion with the simplest 

case (Hi + iH2) / V^, where Hi and H2 are random hermitean Gaussians. The 
essential difference between this case and the case discussed above is that the 
"quark" propagators do not commute, and care must be taken with respect 
to their ordering. 
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The sectors of the wheel (kernel, denoted by F) correspond to the "two-quark" 
irreducible scattering amplitude. In our case it is not 1 (like in the hermitean 
Gaussian case) but a tensor product in "isospin" space 



p/2/3;92S'3 




72/3 




9293 




1 





. (5.6) 



9293 



The two terms correspond to the "gluon" amplitudes of Fig. 3.1b and 3.1c, 
respectively. This particular form becomes transparent in Fig. 5.2. 



91 -(Ql^ 




■h 



■^93 



Fig. 5.2. Two-point kernel with f,g = q,q. 



The factors 
1 

zw 

get transformed (after "dressing") into 



(5.7) 



g{z) g^{w) . 



(5.8) 



Here Q{z) and Q{w) are the resolvents discussed in section 4.1 with r — 0. 
Transposition keeps track of the "flow" of indices. 



Another generalization deals with the treatment of the derivatives. Now the 
logarithm is not a holomorphic functions so we must treat dz as 



dz = ^(9a; - idy) 



(5.9) 



where z — x + iy and 



(5.10) 



where w — u + iv. Note that in order to obtain G{z, z) we cannot write the 
derivatives as d^d^ but rather 8^0^; and set W = ^ at the end of the calculation. 



The choice of "isospin" in the lower fermion loop is done by choosing the 
appropriate derivative d^j for the "quark" and d^j for the "conjugate-quark" . As 
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in [11], two distinct contributions may be identified. First, when the derivatives 
(bullets) act on the exposed or hidden F hues within the same sector of the 
wheel, as shown in Fig. 5.3. Second, when the derivatives act on different 
sectors as shown in Fig. 5.4. 



did2 



Fig. 5.3. Derivatives acting within the same sector of F in the wheel. 



13 + II 



= . . . 9i g^r) . . . 92 g^r) 

Fig. 5.4. Derivatives (bullets) acting on two different sectors of F. 
Finally, we have 



as an apparent generalization of the Brezin-Hikami-Zee formalism. The de- 
nominators in the above formula come from the fact that the dotted lines in 
Figs. 5.3, 5.4 could represent 0, 1, 2, ..oo sectors F of the wheel (see Fig. 5.1b) 
and have to be resummed. The last result may be rewritten as 



N^Gc{z, w) = -a,a^tr,,-iog (i-gi^ g^r) . 



(5.12) 



Here the logarithm is understood as a power series expansion. The operator 
^1 ® ^J'F is a tensor product of 2 x 2 matrices. One can rewrite it as an 
ordinary 4x4 matrix. The trace tr^g becomes an ordinary matrix trace. This 
allows the following reformulation 



N^Gc{z, w) = -dA logdet(l - ® g^F) 



(5.13) 



Using the exphcit expressions for the Ginibre-Girko resolvents (4.1) 



24 




1, the operator Qi (g) can be evaluated to 




In a 4 by 4 matrix notation Qi (g) is 

/ 



01 



^5f^ ^f^w 

ZW QzQw 

g^Qw ZW 
\^ g^w zgy, y 



The determinant oll — Qi® QjF gives — w]^, so 



(5.14) 
(5.15) 



(5.16) 



(5.17) 



N''G,,{z,w] 



(w — zy 



N''G,g{z,w)^0. 



(5.18) 
(5.19) 



5.3 Elliptic Correlator 



We consider now the general case of subsection 4.1. In this case, there is an 
extra contribution to the Ginibre-Girko diagrams, due to the r terms in the 
"gluonic" propagator. In terms of graphs, another diagram may contribute to 
F with a penalty factor r (see Fig. 3.1e and 3.1h). As a result. 



F = 



+ r 



10^ 




^0 


oo) 

















i 




(5.20) 
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The rest of the calculation follows the Ginibre-Girko case, with the resolvents 
for the elliptic case as in section 4.1 



1-t2 



1-t2 



z-TZ 

\ 9z Z-TZ J 
y Qu, W-TW J 



with gl = — r)^ — t{z + z^ — (1 — r^)^- After defining 

Zr = Z — TZ, Wr = W — TW , WC haVC 

N'^Gc{z,w) = -^Alogdetr 



(5.21) 
(5.22) 

(5.23) 



where JC is given explicitly by 



-TZrQu 



Z-rW-r 



-9z9w 

-QzWrW 



-QzWr -TOzQw 
-QzQw -TQzWr 



Z-tWt 



-TZrQu 



Zt — TZt-Wt J 



(5.24) 



After some remarkable cancelations, the determinant L is again simply wp, 
giving the correlators (5.18-5.19). 

Figure 5.5 compares the analytical results (5.18) (solid lines) for the "quark- 
quark" correlator, to the numerical calculations (dashed lines) using an en- 
semble of 50000 matrices of size 100x100 for r = 0.5 (real/major axis is 1.5, 
imaginary/minor is 0.5). The figures on the left are parameterized hy z = 
w — t/\/2{l-\-i), while those on the right are parameterized by 2; = t/\/2{l-\-i), 
w = it. In Fig. 5.6 we show the analytical result (thick solid line) for the 
"quark-conjugate-quark" correlator (5.19) versus the numerical calculations 
for different matrix sizes = 50, 100, 200. The peak around zero in the nu- 
merical calculations is proportional to A^, while its width decreases with larger 
matrices, approaching zero in the large N limit. Our analytical evaluation of 
the correlators disregard delta-functions. 

Let us now investigate the case of real asymmetric matrices defined by (3.5). 
When the 'twisted' propagator 



(5.25) 
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N G.„(z,w) Girko with 'c=0.5 




50 


-50 
-100 
-150 
-200 

50 


-50 
-100 
-150 
-200 



1 1 1 1 1 1 1 


1 1 1 


r 

1 1 1 


2 -1.5 -1 -0.5 0.5 1 1.5 2 
t 

N^Gg^(z,w) Girko with 1=0.5 


1 1 1 


1 1 1 


' — *^ 

1 1 1 





-2 -1.5 -1 -0.5 0.5 1 1.5 2 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

t t 

Fig. 5.5. Numerical (dashed line) versus analytical (solid line) results for the elliptic 
correlators from GUE (see text). 




Fig. 5.6. Same as Fig. 5.5 but for the "quark-conjugate-quark" correlator divided 
by N. The analytical result (5.19) is zero in this region. 

couples at both ends to the same "quark" (as is the case for the self-energy) it 
gives a subleading contribution in (see Fig. 4.1). But in the calculation of 
the correlator it may couple to two different "quark" loops giving a planar con- 
tribution. Now it joins "quarks" of the same species. Similar arguments apply 
for the 'r'-propagator which joins "quarks" and "conjugate-quarks" . Diagram- 
matically we observe that the 'twisted' propagators and the ones considered 
earlier cannot mix in the planar limit so we must add the two resummed con- 
tributions separately — the one calculated earlier and the twisted one given 
by 

-did2 logdet(l - g^{z) ® gj{w) ■ Ft^isted) (5.26) 
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where T twisted is given by the expression obtained from (5.20) after interchang- 
ing 1 and r: 



twisted — 



+ T 




+ 



+ r 




(5.27) 



Here the terms correspond respectively to the "gluonic" amphtudes in Fig. 3.1a, 
3. Id, 3. If and 3.1g. The determinant is now given by 



const ■ {w — z){w — z) . 



(5.28) 



This leads to correlators of the form: 



N'G,,{z,w) = j-^^, (5.29) 
N'G,^{z,w)^ (5.30) 

[w — zy 

Note that these are different expressions than in the complex case, although 
the Green's functions are identical in both cases. 



5.4 Chiral Correlator I (Holomorphic Region) 



In the nonhermitean (chiral) case the correlator in the outside region is cal- 
culated using the same arguments as above. In the outside, the resolvents are 
holomorphic, hence we may consider the derivatives as acting on analytical 
functions. The amplitude T reads in this case 





^01^ 














r = 


















^ooj 


ah 


,0 0, 


cd 




ab 


,10) 



(5.31) 



where ab are the "quark" indices and cd are the "conjugate-quark" indices. 
Note that the kernel is rewritten in the chiral basis, i.e. the indices a,b,c and 
d run from 1 to 2N. The correlator reads 

N'G,{z,z) = -^9,9,logdet(l - Ggg{z) ^ gT.{z)r) (5.32) 
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From the definition of Q in (4.20) we see, tfiat in tlie liolomorpliic case (S2 = 
S3 = and Si = G) the Qgg component (after performing the trivial tracing 
of the diagonal N x N blocks) is simply 



Gqq = 




(5.33) 



where we used the explicit chiral structure (4.19). Performing the inversion 
we arrive at 



1 

det 



z-G -n ^ 

y ^ z-G J 



det\ ^ z-G 



(5.34) 
(5.35) 



with the determinants given by 

det^{z-Gf + fj,^ , dei^{z-Gf + f/ 



(5.36) 



with G being the solution of Pastur equation (4.28) and G is just the complex 
conjugate of G. Now the operator Qi (g) QjT reads 



(5.37) 



1 


1 z-G^ 










( oV 


\det\^ 








+ 






^0 n ^ 


[z-G J 


^z-Go) 



or in 4 X 4 form 



1 

D 



ij,^ i^-cp ^ 

-fi{z -G)00 n{z- G) 
-fi{z -G) ii{z~ G) 
Iz-GP 



(5.38) 



with 



D = \det\ 



G\ 



(5.39) 
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where holomorphic G is the appropriate branch of (4.28) rewritten in the form 
G[{z - Gf + ^J?]^ z-G. The determinant oil-Qi® g^T gives 



^^^^ (5-40) 



The zero of the determinant in (5.32) occurs for [D — jj^) — \z — G\^, that is 

\z-G\'{l-\G\'')- ^''\G\^ (5.41) 

as quoted in [34]. In the case = (and z = w), the determinant in (5.32) is 
simply (1 — G'^{z)G'^{w)) (chiral) as opposed to (1 — G{z)G{w)) (non-chiral) . 
As a result, lor w — z and // = 0, (5.32) is 

N^G{z, z) = l/{z\z^ - Af) (5.42) 



which coincides with (5.5) in [39]. For small z, the behavior G{z,z) ~ l/z"^ 
reflects the exchange of two "massless" modes for 2; = im ^ 0. 

The numerical results for the chiral correlator versus the present analytical 
results are shown in Fig. 5.7. The dashed line is obtained from a numerical 
simulation of 200 100x100 matrices at 11^ = 2, while the solid line follows 
from (5.32). The dashed region indicates the extent of the nonholomorphic 
region for this choice of //^ = 2., over which the result (5.32) no longer holds. 




0.5 1 1.5 2 25 3 

y 



Fig. 5.7. Numerical (dashed) versus analytical (solid line) results for the chiral cor- 
relator I. The nonholomorphic region is shown dashed. 



5.5 Chiral Correlator II (Non-Holomorphic Region) 



As an extension of this construction we consider the case of a nonhermitean 
deterministic matrix plus a hermitean chiral random matrix. The extension 
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amounts to adding an "isospin" index to the resolvents which are now 4x4 
matrices (2 "isospin"+ 2 chirahty). The original 4x4 matrix Qqq{z) (S)^^(w)r 
is now changed into a 16 x 16 matrix (each factor in the tensor product is now 
4 X 4 so the tensor product has dimension 16). The matrix Q{z) ® Q^{w)T is 
of the form 



Qqq 
Qqq 
\Qqq 



Gqq^ 



Gqf G> 
Gjf Oq 



QQ 



^qq 



(5.43) 



The complexity in this case stems from the fact that in the nonholomor- 
phic region "quarks" may turn to "conjugate-quarks" and vice-versa, with all 
"quark" species interacting with themselves. 




X 



Fig. 5.8. Extent of the non-holomorphic region (shadowed) for the chiral ensemble 
at iJp = 0.4. The thick dashed hnes indicate the paths along which the numerical 
correlator was evaluated. 

The resolvents in (5.43) follow from the solutions presented in section 4.3. 
Defining z — x + iy,w — u + iv, we have 



gqq ^(z + II-f- E4) • 



gqq^{z - H'J - A ^ 



El {z - T,4)iy/i^ - jji 

y-{z - T,A)iy / jj. ^ H El 

E4 {z - J:i)iy / n + n 

-{z -Ei)iy/ij,- n E4 



gqg = R-A-^=R 



1 iy/n 
-iy/n 1 



(5.44) 
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N^G (z,w) at n^=0.4 



N^G,„(z,w)atn^= 




1.5 2 



Fig. 5.9. Numerical (thin and dashed lines) versus analytical (thick solid lines) 
results for the chiral "quark-quark" correlators II. The correlators were evaluated 
along the paths shown in Fig. (5.8): z—w=l-2 (left), 2-3 (right). 




Fig. 5.10. Same as Fig. 5.9 for "quark-conjugate-quark" chiral correlator II. 



where 

(z-J:i)(z-^A+li'^-R^^det= , ^ ^ (5.45) 

r -/^^ 



and 



32 



Using (5.44-5.47), the determinant of (5.43) viewed as a 16 x 16 matrix sim- 
plifies to 



det(l-^(.) ^-(^r) = |.-^n.+^r (^^-(^^-^^)(^;-^^))^-"^^^ 5.48) 

A* 



Wc have evaluated the correlator numerically inside the nonholomorphic re- 
gion. Figure 5.8 shows the island for /x^ = 0.4, and indicates the paths z and 
w taken during the numerical evaluation of the correlator. In Fig. 5.9 the an- 
alytical results (5.43) are compared to the numerical simulations. The thin 
solid fine is obtained from a numerical simulation of 73000 200x200 matrices 
at = 0.4, the dashed line represents 75000 events of 100x100 matrices, 
while the thick solid line is deduced from the analytical prediction (5.43). The 
numerical results converge slowly to the analytical one. In Fig. 5.10, the thick 
solid line on the right corresponds to the analytical prediction (5.43), while the 
thin lines on the left are the numerical simulation using 15000 50x50 (dotted), 
15000 100x100 (dashed) and 10000 200x200 matrices (solid line). The height 
of the numerical peak is proportional to A^, diverging in the large N limit. Its 
width is slowly approaching the analytical one. 



5.6 Scattering Correlator 

Here we consider the correlator for the model discussed in section 3.4, that is 

H - i^VV^ ^H- 5VV^ (5.49) 



where H is random Gaussian Hermitean. The precedent discussion can be 
easily adapted to this case by assessing F the "quark-quark" scattering am- 
plitude. In this case, the latter splits into a sum F^ and F^^\ The first one 
is just 

= Sf,f,5g,g, (5.50) 



where /i, /2 are the ingoing and respectively outgoing "quark isospin" on the 
upper line, while the QiS are the corresponding ones for the lower lines (see 
Fig. 5.11a). The graphs F^^^ are considerably more involved (see Fig. 5.11b) 
due to the fact that the V vertices always appear in pairs. Therefore they 
always contain an exchange of two F propagators between the "quark" lines. 
Note that the IPI "quark" scattering amphtude F^^^ flips "quarks" into 
"conjugate-quarks" and vice-versa on the upper and the lower lines. With 
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/l II /2 fl 1 I /2 

S'l— II— gi—i^A\—g2 



^ 9192 



Fig. 5.11. a) Graph contributing to r^^^ (left) and b) graph contributing to 
(right). 



this in mind, 

T 

hf2;9i92 



^VL92 = m ■ F/./. • F,,,, (5.51) 



where the F's are defined by Fig. 4.4 and the factor m comes from the interior 
loop of Fig. 5.11b. Hence {z — x -\- iy ,w — u -\- iv) 



Fqg = d{l + SG)/den (5.52) 
Fgg = -5{1 - SG) /den (5.53) 
Fgg^-S^Ggg/den (5.54) 



where 



den — (5.55) 



and 



Ggg = JGG-- (5.56) 

V 1 -7y 



In the present model the "quark-quark" correlator is now 

N^Ggg{z,w) = -a,a^iogdet(i - g{z) ® g^{w) ■ [r^ + r^^']) (5.57) 



where the dot stands for ordinary matrix multiplication. The determinant in 
(5.57) may be explicitly evaluated to give 



,2(i--y)(i-vO("' + y'')-:.'y^ 



^- • - - < ^;(i-,,)^(i-;^- • ^'-''^ 



As in the case of the elliptic ensemble, the correlators differ depending on 
whether the matrices F and V stem from GUE or GOE ensembles. In the last 
case, a new set of diagrams have to be added. They are leading in GOE and 
subleading in GUE. We have checked that the additional contribution is the 
same as (5.58), resulting into an extra factor of 2 for the GOE correlators. 



34 



Fig. 5.12. Extent of the nonholomorphic region (shadowed) for the scattering prob- 
lem at 5 = 4, m = 0.25. The thick dashed hnes indicate the paths along which the 
numerical correlator was evaluated (see below). 
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Fig. 5.13. Numerical versus analytical (thick solid lines) for the "quark-quark" cor- 
relator of the random scattering problem. 
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Fig. 5.14. Same as Fig. 5.13 but for the "quark-conjugate-quark" correlation function 
divided by N. The analytical result (5.57) is not distinguishable from zero on the 
plot (see text). 
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Figure 5.13 compares the numerical versus analytical results for random scat- 
tering correlator using GUE ensembles along the paths indicated on Fig. 5.12. 
The thick solid line is the analytical result (5.57), while the thin lines follow 
from numerical simulation of 10000 100x100 (dashed), 3000 200x200 (dotted) 
and 2000 400x400 matrices (sohd) at 7 = 4, m = 0.25 along the trajectory 
z = 0.777t(2 + 0.2865i), w = 0.2865it. Part of the path goes outside the non- 
holomorphic region as shown in Fig. 5.12 and the analytical curve is missing 
for that values. The convergence of the numerical results to the analytical one 
is slow. The numerical results for the "quark-conjugate-quark" correlators are 
shown in Fig. 5.14. The analytical result (5.57) is not distinguishable from 
zero in the region shown. The large central peak is a finite size effect, as it 
starts to shrink for N ~ 400. 



6 Partition Functions 



From a "thermodynamical" point of view, the information carried by the one- 
and two-point functions is sufficient to specify the "thermodynamical" poten- 
tial to order 0{1/N) in the entire 2;-plane modulo isolated singularities, as we 
now discuss. Similar ideas were used either in the context of two-dimensional 
models in matter [40], or as a way to calculate next-to-leading corrections to 
hermitean correlators [3], exploiting the functional formalism discussed in [41]. 




Fig. 6.1. Leading contributions to Z]\f in (6.2). 

Let Zn be the partition function in the presence of an external parameter z. 
In the 1/A^ approximation, the diagrammatic contributions to the partition 
function Zn are shown in Fig. 6.1, corresponding to 



log Zj, = NEo + E, + 0{^) (6.1) 

£■0 is the contribution of the "quark" or "conjugate quark" loop in the planar 
approximation, and E-i is contribution of the "quark-quark" loop, and so on, 
in the same approximation. Throughout this section, we will restrict our at- 
tention to nonhermitean matrices with unitary randomness, in which case the 
non- planar corrections to Eq are of order 1/iV^. ^ Hence, Eq is determined by 
the one-point function and Ei by the two-point functions. 



^ A similar observation was made in [3], while evaluating the correction to the 
resolvent for hermitean GUE ensembles. 
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For z such that (6.1) is real, continuous and nondccrcasing function of the 
extensive parameters [42], log Zj^/N may be identified with the "pressure" 
of the random matrix model. As a result, the isolated singularities in the 
"pressure" are just the "phase" boundaries provided that the expansion is 
uniform. Below we give examples where the "phase" boundary is either mean- 
field-driven or fluctuation-driven. 



6.1 Holornorphic Z 

To illustrate the above arguments, consider the partition function 

Zn = {det{z - M)) = Q ##e-^("-^)^^ (6.2) 

In contrast to the one- and two-point correlators discussed above, the deter- 
minant in (6.2) is not singular in the ^-plane configuration by configuration. 
Hence, (6.2) is a priori holomorphic in z (minus isolated singularities). 

The one- and two-point functions on their holomorphic support may be ob- 
tained from logZjv by differentiation with respect to z. Therefore, from (6.1) 

2 

Eq^ j dz' G{z') + canst (6.3) 

or equivalently 

Eq^zG- jdC z{G) + const (6.4) 

after integrating by part. Note that z{G) = B{G) is just the Blue's func- 
tion [31] of G (where the Blue's function is the functional inverse of the resol- 
vent, i.e. B[G{z)\ — z). The constant in Eq is fixed by the asymptotic behavior 
of (6.2), that is Z^ ~ z^ . The planar contribution to Ei in (6.1) follows from 
the "quark-quark" wheel of Fig. 6.1, that is 

El = -^log det(l -g®g^ -T) (6.5) 

which can be understood by opening one of the lines in the two-point correlator 
(5.12): the factor 1/2 is combinatorial (two identical "quarks") and the minus 
sign corresponds to decreasing the number of fermion loops by one. The deriva- 
tive amputates the left out propagator after "opening" the fermionic loop. We 
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note that the factor —1/2 in (6.5) causes the "quark-quark" contribution to 
be overall "bosonic" ^ . Hence 

Zn = e^^" ■ ({det(l -g®g^ -V)]''^ +0 Q)) (6.6) 

In contrast to the "quark" contribution, the "bosonic" contribution is not 
extensive in N , since in random matrix models is like "color" not "space" . 
The "bosonic" contribution is dwarfed by the "quark" contribution (1 : A^) 
[43] . Both Eq and Ei are simple functions of the resolvent on a specific branch 
as we will show in the examples below. When the "quark-quark" contribution 
(pre-exponent) to (6.6) does not vanish in the 2;-plane, the singularities of 
are those of the "quark" contribution in Eq in large N . Such singularities are 
mean-field-driven. Otherwise, they are fiuctuation-driven. 

We note that the partition function Zjv through (6.6) exhibits an essential 
singularity in 1/N as expected from thermodynamical arguments. Assuming 

that the expansion for \ogZ]^/N is uniform, then logZ^v/A^ follows from (6.6) 
using the holomorphic resolvent G{z) for large z. The small z region follows by 
analytical continuation. However, since G{z) is multi-valued (already the sim- 
ple case of the Ginibre-Girko ensemble yields two branches for the resolvent 
in (4.7)), the analytical continuation is ambiguous. The ambiguity may be re- 
moved by identifying log \ Z]^\/N with some generalized "pressure" and taking 
G{z) so that log \Zm\/N is maximum ^ . As a result, log Z^/N is piece-wise 
analytic in leading order in with "cusps" at 

F(^^)(x, y) = V{}\x, y) - V^\x, y) = , (6.7) 

following the transition from branch i to branch j of G with V^r = log \ Zi^\/N. 

The character of the transition in the 1/N approximation can be highhghted 
by noting that for any finite N , the partition function (6.2) is a complex 
polynomial in z of degree N with random coefficients. In large N , 

1 If 
= — log |Zjv| = - / dvdvQ{v,v)\og\z-v\^ . (6.8) 



To leading order, the distribution of singularities along the "cusps" (6.7) is 



^ A similar result can be obtained using a bozonized version to (6.2) in the Gaussian 
approximation. 

^ This condition is equivalent to the saddle point condition. It does not necessarily 
fulfill the conditions discussed in [42] . 
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Q{z,z) = ^\d,F\' S{F{z,z)) (6.9) 

which is normahzed to 1 in the z-plane. Redefining the density of singularities 
by unit length along the curve F{z, z) — 0, we may rewrite ( 6.9) as 



e{z,z)\^^, = ^\d,F\ ^ - G'(^-)| . (6.10) 

For g ^ 0, the integrand in (6.8) is singular at z — v which results into different 
forms for V^, hence a cusp. For g = 0, that is dzF = 0, is differentiable. For 
physical Vjy (real and monotonically increasing), the points g = are multi- 
critical points. At these points all n-point (n > 2) functions diverge. This 
observation also holds for Ising models with complex external parameters [44] . 
Assuming macroscopic universality [10] for all n-points {n > 2), we conclude 
that dzF = means a branching point for the resolvents, hence dzG = oo 
or B'{G) = [31]. For hermitean matrices, these conditions coincide with the 
end-points of the eigenvalue distributions [31,32]. 



• Example 

To illustrate the above concepts, consider the Ginibre-Girko ensemble of sec- 
tion 4.1. The resolvent in the holomorphic region satisfies (4.7), so 



z^tG + ^. (6.11) 

The integration (6.4) in Eq is straightforward, and after fixing the asymptotic 
behavior we obtain 



= G-^e5^«^ ((1 - G'{z)r)-h + (^)) ■ (6.12) 

Here G is the solution of (6.11). The pre-exponent in (6.12) follows from (6.6) 
with the matrix Q replaced by G and F = r, as seen in the "quark-quark" 
component of the vertex matrix in (4.1). Using (6.11) we observe that the 
pre-exponent diverges at two points in the z-plane, z"^ — At. At these points 
there is a "phase" change as we now show. 

Given (6.12), the generalized "pressure" in leading order is 



V^^-hog{G^G^) + ^{Gl + Gl) + o(^^] . (6.13) 
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V± define two intersecting surfaces valued in the 2;-plane, for two branches G± 
of the solutions to (6.11). The parametric equation for the intersecting curve 
is 



F{z,z) ^V+-V_^0 (6.14) 

As indicated above, Vj^ is piece-wise differentiable. Note that F = on the 
cut along the real axis, — 2i/t < z < +2^/t, and from (6.10) the density of 
singularities per unit length is 



1 Vir" 



z 



2 



TT ZT 



Along F, the density of singularities is semi-circle. The density (6.15) vanishes 
at the end-points z = ±2^^- This is easily seen to be the same as d^G = oo, 
or dB{G) / dO = with B{G) = tG + 1/G. As noted above the term in 
bracket in Eq. (6.12) vanishes at these points, with a diverging "quark-quark" 
contribution. The transition is fluctuation-driven. These points may signal 
the onset of scaling regions with possible universal microscopic behavior for 
nonhermitean random matrix models. This issue will be pursued elsewhere. 
At these points the 1/A^ expansion we have used breaks down. 

Note that the case of the circular ensemble with r = is very specific. Prom 

our analysis, we see that the support of F shrinks to a single point z = 0. In 
fact Zjsf — z^ to all orders in 1/N. This can be understood by noting that 
since in 

Zjv = (det(z - M)) (6.16) 



the determinant is just 

det(^ -M)= e,^,..,,{z - M)l{z - M)l ...{z- M)^ (6.17) 



a// expectation values vanish (exactly to all orders in 1/A^) for random gaussian 
complex matrices, i.e. 

(M?M^...M^) = 0. (6.18) 



Hence our claim. 

The present construction for the partition function for the Ginibrc-Girko en- 
semble extends to the other cases of nonhermitean ensembles with unitary 
randomness. We briefly mention the case of chiral random ensembles plus 
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deterministic nonhermitean matrix of (off-diagonal) strengtli The relevant 
one- and two-point correlators in the holomorpliic and nonholomorpliic do- 
mains were explicitly constructed in sections 4.3, 5.4 and 5.5. Using these 
results and in analogy with (6.6), elementary integration for this case leads to 

Z^(z,C) = e^^°-({i^-^[(i^ + CY-(^-G)l}~^+oQ)) , (6.19) 

where D = {z - Gf/G^, and 

E,{z,C) = G' + log^ (6.20) 

with the appropriate branch of holomorphic G solution to (4.28), with fJ^ = (. 
Note that for 2; = and G^ = —1 — C^, the pre-exponent in (6.19) diverges. 
It also diverges at z — which are the zeros of (6.10) for the present case 
(two zeros for small ( and four zeros for large (). Again, at these points, 
the 1/N expansion breaks down marking the onset of scaling regions and the 
possibility of microscopic universality. The z = divergence is just the notori- 
ous "Goldstone" mode in chiral models, illustrating the noncommutativity of 
— > 00 and 2; — > 0. The rest of the arguments follow easily from the preceding 
example, in agreement with the "thermodynamics" discussed in [43]. 

• Replicas 

We note that in the circular case, the partition function with Nf "replicas" 

^iv,iv, = (det^/(^-M)) (6.21) 

yields simply 

Zn,n, = z'^^'f (6.22) 

For general Ai and Nf noninteger the determinant in (6.21) is multivalued, and 
a determination is in general needed. Here, it will be assumed by interpolating 
the integer result [43] ^ . While the limits A" — > 00 and Nf ^ are compatible 
for In Zj^^fqjN through the diagrammatic analysis, it is not necessarily the case 
for the connected functions (^-derivatives) since singularities for 2; ~ are 
generated. 

^ Another way is to use the ^-function in combination with the heat-kernel con- 
struction, through det^t A = e'^f^^W and ^a{s) = Tr(l/^*). We have checked 
that the leading term is in agreement with the interpolation for smooth operators. 
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In the expansion, the "quark" insertions are suppressed by powers of 
Nf/N for singularity- free operators. In general, this implies that random ma- 
trix models are self-quenched [43] , and holomorphy in the external parameter 
z is preserved as — > oo (modulo cusps). In the presence of singularities, 
holomorphic symmetry is usually upset as A?" — > oo for a certain range of Nf 
and z A4, and arguments similar to those of section 3 should be used. This 
confirms the shortcomings of (6.21) to yield the quenched one-point functions 
[23], and two-point functions [34] for small z. 



6.2 Nonholomorphic Z 



The above analysis for the holomorphic thermodynamical potential may also 
be extended to nonholomorphic partition functions, again for Gaussian unitary 
randomness. For example, 

Zn[z, z] = (det \z -M\^) = l^j ####e-^(^--^)'^-^(^-^^)'^^ (6.23) 



which is the random matrix analogue of the partition function discussed in [45] 
for unquenched QCD. In these partition functions the phase of the {Nf — 2) 
determinant is set to one. Again, the one- and two-point functions may be 
obtained from log by differentiation with respect to z and z (modulo delta 
functions) . 

Prom (6.1) we have 



G{zrz) 



dEo{z,z) 
dz 



G{z,z)^ 



dEo{z,z) 
dz 



(6.24) 



These expressions may be integrated and the integration constant is set by 
the matching condition 



Eq ~ log 1^1 



for 2; — >• oo . 



(6.25) 



For the Ginibre-Girko ensemble the nonholomorphic resolvent (4.2) yields 



Eo{z,z) = 



1-t2 



+ Cnh 



(6.26) 



while the holomorphic resolvent (4.7) gives 



z r 1 

Eo{z, z) ^ — (z^ Vz'^ - 4r) ± log -{z + V z^ - 4t) 
4t ^ ^ ^ 



+ h.c. + Ch (6.27) 
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with the hmit 



hm Eq{z, z) ^1 + log \zf + Ch . 



(6.28) 



The constant c/j is — 1, and Cnh is determined from the continuity condition of 
Eo at the boundary (4.9). At r = 



Eq 



log 

< 


\z\^ • 


\z\ 






\z? 


- 1 : 


\z\ 



(6.29) 



The "phase" change at = 1 is fluctuation-driven. 



The expression for the correlator in the holomorphic phase follows by differ- 
entiating the term E\^ that is 



d^d^Ex = (tr 



1 



z-M 



(6.30) 



Hence 



Ei^ - \ndet{l-g®g • T) + $ 



(6.31) 



where $ is a harmonic function. To determine $ we note that 

d,E,^Gi{z), 



(6.32) 



Gi{z) is the correction to the one-point function in the planar approx- 
imation. It follows from the insertion of "quark" or "conjugate-quark" loops 
in G{z). Specifically 



which can be rewritten as 



—T 1 

•r) + -lndet(l-^®^^ 



|lndet(l - 
A similar expression follows from G. Hence 



D} 



(6.34) 



El = iogdet(i -g^g^ -r)-^ iogdet(i -g^g^ -r) 
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logdet(l-e? (8) • r) - ^ logdet(l-^? -F). (6.35) 
As a result, the partition function (6.23) reads 

Z^[^,^]=e^^° • (j^det{l-g 0^ ■T)\det{l-g 0^ ■T)\y + o(^^J^ .(6.36) 

Note that the contributions from the wheel-diagrams are of the form 1 / V det 
and hence "bosonic" in character. In (6.35), there are two contributions from 
the "quark-conjugate-quark" wheels (upper line in (6.35)), one contribution 
from the "quark- quark" wheel and one contribution from the "conjugate- 
quark-conjugate-quark" wheel, thereby explaining all contributions in (6.36) 
at next to leading order. Like in the holomorphic case, again similar results 
can be obtained using a bozonized version to (6.23) in the Gaussian approxi- 
mation. 

Again, the partition function Zj^ has an essential singularity in 1/N, but 
log Zn/N does not. For any finite N, the latter diverges for 

det(l - ^ ® . r) = (6.37) 

The divergence condition from the second pre-exponent in (6.36), i.e. 

det(l -g^g^ ■T) = (6.38) 

is not essential, since the line of singularities following from (6.37) defines the 
boundary of a surface that in general contains the cusps at dzG = oo. In 
particular, the set of discrete points (6.38). In this case, the transition in the 
2;-plane is fluctuation-driven. For the Ginibre-Girko example, the nonholomor- 
phic partition function diverges according to (6.37), that is for 

1 - = (6.39) 

Indeed in the holomorphic region, we could replace g by G, the modulus 
follows from the "conjugate-quark" resolvent, and the vertex corresponds to 
the "quark-conjugate-quark" element of the kernel (4.1) set to 1. The line of 
singularities (6.39) reproduces in this case the ellipse (4.9). The ellipse includes 
the points of the "phase" change, i.e. the focal points = 4r, connected by 
the interval (6.14), i.e. F = 0. We recall that these focal points were obtained 
as a divergence condition for (6.12), a special case of the general condition 
(6.38). 

The same results were also discussed in [34] for the case of a chiral random 
matrix model with a chemical potential /i — in the quenched approximation. 
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In particular, it was shown that the condition (6.37) exactly reproduces the 
islands of "mixed-condensate" , obtained using the replica methods [23]. In this 
case the condition (6.37) reads [34] 

D-''[{D - C')' -\z- G\^] = (6.40) 

with D = \{z — G)/G\'^, therefore exactly the condition (5.41). Comparing to 
the "quark-quark" correlator in (6.19), we notice the appearance of the mod- 
ulus and the flip in the sign of ( due to the fact that now we measure the 
correlation between the "quark" and the "conjugate-quark". Equation (6.40) 
defines the boundaries of the fore mentioned "islands" . Inside (6.40), the resol- 
vents in (6.36) are nonholomorphic, while outside (6.40) they are holomorphic. 
The holomorphic and nonholomorphic resolvents match on (6.40). 



7 Conclusions 

We have shown that chiral and hermitean random matrix models are amenable 
to diagrammatic techniques in the context of the 1/A^ approximation. Since 
the matrices are hermitean, the eigenvalue distributions and their macroscopic 
correlations follow from one- and two-point functions that are analytic in 
the complex z-plane modulo isolated singularities. Although our results were 
carried using "quenched" measures, we expect them to hold for "unquenched" 
ones as well since the "quark" effects are 1/A^ suppressed [43]. 

We have shown how to extend the diagrammatic arguments of Brezin, Hikami 
and Zee to a large variety of nonhermitean random matrix models. For non- 
hermitean matrices, the one- and two-point correlation functions are analytic 
in the complex z-plane modulo surface singularities. These singularities origi- 
nate from an accumulation of eigenvalues in the large N limit, and cause the 
one- and two-point functions to be nonholomorphic in finite domains. The 
breakdown condition is set by the divergence of the correlations between pairs 
of eigenvalues, conjugate of each other and a distance N apart in the spectrum 
[34]. We have exphcitly constructed the one- and two-point functions in the 
planar approximation and shown that the analytical results are reproduced by 
direct numerical calculations using large samples of random and nonhermitean 
matrices. 

We have explicitly shown how the partition functions of nonhermitean ma- 
trix models with unitary randomness can be related to the one- and two-point 
functions for holomorphic as well as nonholomorphic kernels to order 0{1/N), 
thereby providing a simple physical interpretation to several preceding results. 
The partition functions to the order calculated are solely a function of the re- 
solvent G. We expect this property to hold to higher orders, and even suspect 
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a generalized form of macroscopic universality when extended to non-Gaiissian 
weights [10]. Generic equations for the singularities of the partition functions 
were found both for the holomorphic as well as the nonholomorphic cases 
to order 0{1/N). These singularities are related to structural changes in the 
random matrix model, such as changes in the eigenvalue distributions and 
correlations in models with [23.34] or without dissipation [32,43]. Our con- 
struction may be useful for investigating a number of related issues in random 
matrix and polynomial models [46]. 



Acknowledgements 

This work was partially supported by the US DOE grant DE-FG-88ER40388, 
by the Polish Government Project (KBN) grant 2P03004412 and by the Hun- 
garian Research Foundation OTKA grants T022931 and F019689. The au- 
thors would like to thank Jurek Jurkiewicz for discussions. I.Z. thanks Robert 
Schrock for a discussion. G.P. thanks the Soros Foundation for financial sup- 
port and the Nuclear Theory Group at Stony Brook for hospitality. R.A.J, 
would hke to thank the Nuclear Theory Group at Stony Brook and GSI, 
where part of this work was done. 



Note Added 

After submission of this work, some of our results were confirmed by other 
authors. Our diagrammatic approach and part of our results have been con- 
firmed in two interesting papers by Feinberg and Zee, cond-mat/9701387 and 
cond-mat/9704191. The analytical results of section 6 as illustrated in our 
second example and the discussion on the replicas, account for the nature 
and location of all singularities and their distribution as they appeared in 
an extensive numerical analysis (up to 500 digits accuracy) by Halasz et al., 
hep-lat/9611008 and hep-lat/9703006. 
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